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Abstract — We review the main aspects of the fragmentation of bottom quarks in top quark decay. 
The NLO b-quark energy spectrum presents large mass logarithms ln(mf/m|), which can be re- 
summed by the use of the approach of perturbative fragmentation functions. Large soft contributions 
in both coefficient function and initial condition of the perturbative fragmentation function have been 
resummed as well. Results on the energy distribution of b quarks and b-flavoured hadrons are finally 
presented in both x and moment spaces. 

1 Introduction 

A reliable understanding of bottom quark fragmentation in top quark decay (t — > bW) will be fundamental to 
accurately measure the top properties, such as its mass m t , at present and future high-energy colliders. In fact, 
the uncertainty on bottom quark fragmentation is one of the sources of systematic error on m t at the Tevatron 
accelerator 1 1 1 and will play a crucial role in the reconstruction of m t from final states with leptons and J/ipat the 
LHC0. 

In this paper we investigate bottom fragmentation in top decay within the framework of perturbative fragmentation 
functions [3 1. We shall resum collinear logarithms ~ ln(m| / mT) and soft terms that appear in the next-to-leading 
order (NLO) b-quark energy distribution. We shall present results on the b-quark energy spectrum in top quark 
decay and investigate the impact of collinear and soft resummation. Hadron-level results on b-flavoured hadrons 
will be shown in xb and moment spaces. 



2 Collinear and soft resummation 

InRef. |4 1 NLO corrections to top decay t{pt) — > b(pb)W(pw)(g(Pg)) have been computed for a massive b quark, 
and the differential width dT/dxb, with x b being the normalized b-quark energy fraction in the top rest frame, has 
been calculated. The differential rate obtained in |4 1 exhibits large mass logarithms ~ as ln(mj /m?) that need to 
be resummed in order to improve the prediction. 

Such contributions can be resummed by using the perturbative fragmentation approach which, up to power cor- 
rections, factorizes the rate of heavy-quark production into the convolution of a coefficient function, describing the 
emission of a massless parton, and a perturbative fragmentation function D(/j,p,m), where [ip is the factorization 
scale. In the MS factorization scheme we have: 
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In Eq. Q To is the Born width of the process t — > bW. The 0(a>s) top decay coefficient function has been 
computed in |4|. 



The perturbative fragmentation function expresses the transition of the massless parton into the massive quark, and 
its value at any scale /i p can be obtained by solving the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
evolution equations |5 6 \ once an initial condition at a scale flop is given. 

In the NLO expression for D([iof, m), which was argued to be process independent, was given. The process 
independence has been lately established in a more general way in Ref. Q. 

The initial condition of the perturbative fragmentation function reads [3 j: 



D p {xb,HoF,m b ) = 6(1 - Xb) + 



2tt 



ln^of -21n(l-a:&)-l 
1 - xb V m b 



(2) 



As discussed in |4|, solving the DGLAP equations for the evolution flop — > /xp, with a NLO kernel, allows one to 
resum leading logarithms (LL) ~ ctg \n n (fi F /u,Q F ) and next-to-leading logarithms (NLL) ~ Og In™ -1 (wp / fiQ F ) 
(collinear resummation). If we set u,qf ~ nib and up ~ m t , we resum large logarithms ~ ln(m 2 /m 2 ), which are 
indeed the terms appearing in the massive, unevolved, NLO dT/dxb. 

Moreover, both the MS coefficient function |4| and the initial condition of the perturbative fragmentation function 
present, at 0(as), terms that behave like 1/(1 — Xb)+ or [ln(l — Xb)/(1 — £&)]+, which become large for 
Xb — > 1, i.e. for soft-gluon radiation. In Mellin moment space, such contributions correspond to behaviours 
~ In N and ~ In 2 N respectively. 

Soft contributions in the perturbative fragmentation function are process-independent and have been resummed in 
f7) with NLL accuracy. Soft terms in the coefficient function are instead process-dependent. Resummation of LL 
~ ln n+1 N and NLL ~ In™ N contributions to the top-decay coefficient function has been performed in 1 8 1 
and we do not report here the formulae for the sake of brevity. 



3 Parton-level results 

We would like to present results on the b-quark energy distribution in top decay and investigate the effect of 
collinear and soft resummation. In Fig.^we show the b-quark energy spectrum. The NLO calculation lies below 
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Figure 1: b-quark energy distribution in top decay, according to the unresummed fixed-order calculation (dotted 
line), and after inclusion of collinear resummation (dashed) and of both collinear and soft resummations (solid). 
We have set nit — 175 GeV, m& = 5 GeV, mw = 80 GeV, Aqcd = 200 MeV. In the inset figure, we show the 
same curves on a logarithmic scale, for Xb > 0.8. 

the two resummed predictions and is divergent as Xb — » 1. After the resummation of collinear terms ~ ln(m 2 /m 2 ) 
the distribution exhibits a sharp peak at Xb close to 1 . Finally, the inclusion of soft-gluon resummation smoothens 
out the distribution, which exhibits the so-called Sudakov peak. 

As discussed in 1 8 1, the implementation of collinear and soft resummation leads to a milder dependence of observ- 
ables on the factorization and renormalization scales entering the calculation, which corresponds to a reduction 
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of the theoretical uncertainty. As an example, Fig. [2] shows the results on the dependence of the x b spectrum on 
the factorization scale /ip in Eq. Q, which is taken equal to m t /2, m t and 2m t , and the effect of soft resumma- 
tion. We note that while the unresummed prediction still exhibits a dependence on the value chosen for fj,p, the 
implementation of soft resummation yields three almost undistinguishable distributions. 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 2: b-quark energy spectrum for different values of the factorization scale [If, with (solid) and without 
(dashes) NLL soft-gluon resummation. 



4 Hadron-level results 

We would like to make predictions for the spectrum of b-flavoured hadrons in top decay. We write the normalized 
rate for the production of B-hadrons B as a convolution of the rate for the production of b quarks and a non- 
perturbative fragmentation function D np (x): 

ldT B 1 f 1 dzdT» s n np(XB\ ... 

^ — {xB:mt,m b j = - —(z,m t ,m b )D p [ — I, (3) 

1 dxB 1 J XB z dz \ z J 

where xb is the B normalized energy fraction. The parton-level rate dT b /dz can be computed following the 
method which has been discussed in the previous section. 

As for the non-perturbative fragmentation function, one can use some phenomenological models with tunable 
parameters, which are to be fitted to experimental data. We consider a power law with two parameters: 

D np (x;a,P)= 1 (l-x) a xP, (4) 

B(p + 1, a + 1) 



the model of Kartvelishvili et al. (9) 



D np (x;S) = (l + tf)(2 + <*)(! -x)x 5 (5) 



and the non-perturbative fragmentation function of Peterson et al. HI 01 : 

D ^ {X] E) = x[l-l/x A e/(l-x)r (6) 

In Eq. 0, B(x,y) is the Euler beta function; in (|6j A is a normalization constant. We tune such models to e + e~ 
data from the ALEPH (TJ1 and SLD (13 Collaboration. The ALEPH data refer to b-flavoured mesons, the SLD 
data to baryons and mesons. When we do the fits, we must describe the e + e~ — > bb process within the same 
framework as we did for top decay, i.e. we use the perturbative fragmentation method, NLL DGLAP evolution and 
NLL soft resummation. As in |4j[8), we shall consider xb values within the range 0.18 < xb < 0.94. 

In Table^we show the results of our fits, along with the corresponding values of x 2 per degree of freedom. One 
can see that the power law with two parameters ® and the Kartvelishvili model Q fit both ALEPH and SLD 
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Table 1 : Results of fits of hadronization models to ALEPH and SLD data on b-fiavoured hadron production in 
e + e~ annihilation. 





ALEPH 


SLD 


a 


0.51 ±0.15 


2.04 ±0.38 


P 


13.35 ±1.46 


25.18 ±3.27 


X>,/3)/dof 


2.56/14 


11.50/16 


5 


17.76 ±0.62 


16.59 ±0.49 


X 2 (S)/dof 


10.54/15 


22.19/17 


e 


(1.77 ±0.16) x 10~ ;5 


(1.61 ±0.14) x 10~ 3 


X*(e)/dof 


29.83/15 


158.15/17 



data rather well, while the Peterson fragmentation function is marginally consistent with ALEPH and unable to 
reproduce the SLD data. Moreover, the values of the best-fit parameters 5 and e, fitted to ALEPH and SLD, are in 
agreement within two standard deviations. 

In Fig. |3] we show our prediction for the B-hadron spectrum in top decay, using all three hadronization models 
fitted to the ALEPH data. In order to account for the uncertainties on the best-fit parameters, for each model we 
plot a band corresponding to a prediction at one-standard-deviation confidence level. From Fig.|3]we learn that the 
predictions based on the models (0J and are consistent, while the Peterson model yields a distribution that lies 
quite far from the other two and is peaked at larger values of xb- 

In Fig. 0] we plot the xb spectra yielded by models (0} and (0, but fitted to SLD. Such distributions statistically 
agree at the confidence level of two standard deviations. 

In Fig. |5] we compare the predictions obtained using the power law with two parameters, but fitted to ALEPH and 
SLD data. We observe that the spectra are distinguishable; this difference may be related to the different hadron 
types that the two experiments have reconstructed. 




Figure 3: B-hadron energy spectrum in top decay, according to the power law (solid line), the Kartvelishvili 
(dashed) and the Peterson model (dotted), fitted to the e + e~ — > bb data from ALEPH. The plotted curves are the 
edges of bands at one-standard-deviation confidence level. 

We finally wish to present results on the moments of the B-hadron spectrum in moment space Tjy. Such moments 
can be written as the product of a perturbative and a non-perturbative contribution = T^D 7 ^ . The advantage of 
working in moment space is that one can extract D 7 ^ from e + e~ data without relying on any specific hadronization 
model. 

Predictions for the moments Tfj of B-meson spectra in top decay are given in Table|2 where data from the DELPHI 
Collaboration lf]~3l are used to obtain the non-perturbative information D 7 ^. Two sets of perturbative results ([A] 
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Figure 4: As in Pig. O but fitting the hadronization models in Eqs. © and (jSJi to the SLD data. 




Figure 5: B spectrum in top decay according to the power law fitted to ALEPH (solid) and SLD (dashes) data. 

and [B]) are shown, the first using Aqcd = 0.226 GeV and mi, — 4.75 GeV, the second Aqcd = 0.2 GeV and 
mb = 5 GeV, the default values of this analysis. As expected, the perturbative calculations and the corresponding 
non-perturbative components differ at the level of few per cent, according to whether one uses set [A] or [B]. 
However, the final hadron-level predictions for the physical results differ only at the level of per mille. 

5 Conclusions 

We have considered bottom quark fragmentation in top quark decay t — > bW. We have pointed out that the 
fixed-order result on the b-quark energy spectrum exhibits large mass logarithms ~ ln(mf/m^), which can be 
resummed to NLL accuracy using the approach of perturbative fragmentation functions and DGLAP evolution 
equations. Moreover, NLL soft contributions to the coefficient function and to the initial condition of the perturba- 
tive fragmentation function have been resummed as well. 

We have presented results on the b-quark spectrum in top decay, which displays a remarkable impact of the inclu- 
sion of soft and collinear resummation. In particular, the distributions exhibit very little dependence on factoriza- 
tion and renormalization scales. 

Predictions on b-flavoured hadron energy distributions in top decay have been obtained using ALEPH and SLD 
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data to parametrize some hadronization models in xb space and DELPHI data to get non-perturbative information 
in moment space. 

The considered approach can now be applied to study several observables, which are relevant to top quark phe- 
nomenology at the Tevatron and ultimately at the LHC and compare the obtained results with the ones given by 
Monte Carlo event generators. 

Table 2: DELPHI data for the moments a^, the resummed e+e" perturbative calculations for cr b N Q, and the 
extracted non-perturbative contribution D 7 ^ , Using the perturbative results T b N , a prediction for the physical 
observable moments Y% is given. Set [A]: A QC d = 0.226 GeV and m b = 4.75 GeV, set [B]: A QC d = 0.2 GeV 
and m;, = 5 GeV. 





(x) (x 2 ) (x d ) (x 4 ) 


e + e data 


0.7153±0.0052 0.5401 ±0.0064 0.4236±0.0065 0.3406±0.0064 


e+e" NLL a b N [A] 
e+e" NLL a b N [B] 


0.7666 0.6239 0.5246 0.4502 
0.7801 0.6436 0.5479 0.4755 


D7 [A] 
D% [B] 


0.9331 0.8657 0.8075 0.7566 
0.9169 0.8392 0.7731 0.7163 


t-decay NLL T b N [A] 
t-decay NLL T b N [B] 


0.7750 0.6417 0.5498 0.4807 
0.7884 0.6617 0.5737 0.5072 


t-decay T% [A] 
t-decay r$ [B] 


0.7231 0.5555 0.4440 0.3637 
0.7228 0.5553 0.4435 0.3633 
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